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Abstract. We constrain the effective theory of one-body dark matter-nucleon interactions 
using neutrino telescope observations. We derive exclusion limits on the 28 coupling constants 
of the theory, exploring interaction operators previously considered in dark matter direct 
detection only, and using new nuclear response functions recently derived through nuclear 
structure calculations. We determine for what interactions neutrino telescopes are superior 
to current direct detection experiments, and show that Hydrogen is not the most important 
element in the exclusion limit calculation for the majority of the spin-dependent operators. 
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1 Introduction 

The clustering of dark matter in large spheroidal halos hosting galaxies and other visible 
astrophysical structures has inspired complementary approaches to reveal the nature of dark 
matter [1]. The indirect detection of energetic neutrinos from dark matter annihilation in 
the Sun is a prime example of strategies designed for discovering dark matter in astroparticle 
physics [2]. 

Dark matter particles from the Milky Way dark matter halo might be captured and an¬ 
nihilate in the Sun, if they loose energy while crossing and scattering in the solar medium [3] . 
Energetic neutrinos from dark matter annihilation in the Sun are expected in a variety of 
models for weakly interacting dark matter [4-7], and could potentially be observed on Earth 
at neutrino telescopes [8, 9]. 

The predicted flux of dark matter-induced neutrinos from the Sun depends on how 
dark matter is distributed in the solar neighborhood, on the primary dark matter annihi¬ 
lation channel and total annihilation rate, and on the cross-section for dark matter-nucleus 
scattering in the Sun [10, 11]. 

In this work we focus on the cross-section dependence of the expected dark matter- 
induced neutrino flux, in a comprehensive analysis of a variety of dark matter-nucleon inter¬ 
actions. 

We calculate the flux of neutrino-induced upward muons from dark matter annihilation 
in the Sun at neutrino telescopes in the general effective theory of one-body dark matter- 
nucleon interactions. We then compare our predictions with current neutrino telescope ob¬ 
servations, deriving exclusion limits on the isoscalar and isovector coupling constants of the 
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theory. We also compare our limits from neutrino telescopes with current limits from dark 
matter direct detection experiments, establishing when the former are superior to the latter. 
For the elements in the Sun, we consider new nuclear response functions which recently ap¬ 
peared in the literature [12], and analyze the contribution of different isotopes to the final 
exclusion limits in detail. 

Compared to previous analyses, here we explore a larger set of dark matter-nucleon 
interactions, including several momentum or velocity dependent operators previously con¬ 
sidered in dark matter direct detection only. Most importantly, for the dark matter-nucleon 
interactions studied here, and for the 16 most abundant isotopes in the Sun, we consider 
new nuclear response functions recently computed through numerical nuclear structure cal¬ 
culations in [12]. The use of improved nuclear response functions strengthen our conclusions, 
which do not rely on the approximate exponential form factors often used in the literature 
in the analysis of neutrino telescope data. 

Effective theories for dark matter-nucleon interactions were explored in the context of 
dark matter direct detection in [13-32]. Neutrino telescope observations were interpreted in 
terms of spin and momentum/velocity dependent dark matter-nucleon interactions in [33], 
and more recently in [34]. Both analyses focus on a subset of operators studied here, and 
consider either dark matter scattering from Hydrogen only, or simple Helm form factors. 
The effective field theory interpretation of searches for dark matter annihilation in the Sun 
in [35] considers 4 operators, and constant spin-dependent cross-sections. Recently, advances 
in the theory of dark matter heat transport in the Sun within models with momentum or 
velocity dependent scattering cross-sections, and in the context of Helioseismology were done 
in [36-38]. 

The paper is organized as follows. In Sec. 2 we review the general effective theory of 
one-body dark matter-nucleon interactions. We compare the predictions of this theory to 
physical observables introduced in Sec. 3 as explained in Sec. 4, where we also interpret our 
findings. We conclude in Sec. 5. Finally, we describe our LUX and COUPP analysis in 
Appendix A, and list key equations in Appendix B. 

2 Theoretical framework 

In this section we review the general effective theory of one-body dark matter-nucleon inter¬ 
actions mediated by a heavy spin-1 or spin-0 particle. We compare the predictions of this 
theory with current neutrino telescope observations in Sec. 3. 

2.1 Effective theory of dark matter-nucleon interactions 

We start with an analysis of the amplitude A4 for non-relativistic elastic scattering of a dark 
matter particle of mass and initial (final) momentum p (p') from a single nucleon of 
mass mi\f and initial (final) momentum k (k'). The amplitude A4 is restricted by Galilean 
invariance and momentum conservation. Momentum conservation implies that only three 
out of the four momenta p, p/ k and k' are independent. A possible choice of independent 
momenta is p, k and q = k — k', where q is the momentum transfer. Galilean invariance 
constrains A4 to be a function of the initial relative velocity v = pjm^ — k/m 7 v, rather than 
of the momenta p and k separately. Therefore, the scattering amplitude A4 can in general 
be written as a function of q and v, and of the nucleon and dark matter particle spins, Sv 
and S^, respectively. 
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Table 1. Complete set of non-relativistic quantum mechanical operators that are at most linear in 
the transverse relative velocity operator V'*-, and in nucleon and dark matter particle spin operators, 
Sat and S^, respectively. Introducing the nucleon mass, mjv, in the equations above all operators 
have the same mass dimension. 


The non-relativistic quantum mechanical Hamiltonian density that underlies the ampli¬ 
tude M must have the form [16] 

= E (2.1) 

T=0,1 k 

where r is the dark matter-nucleon relative distance, = 1 is the identity in isospin space, 
and = r 3 is the third Pauli matrix. The quantum mechanical operators Ok in Eq. (2.1) are 
constructed from the momentum transfer operator q, the transverse relative velocity operator 
v-*-, and the nucleon and dark matter particle spin operators Sjv and respectively. If we 
demand that the operators Ok are at most linear in v-*-, Sjv and S^, then only 14 independent 
operators appear in the linear combination (2.1)^. They are listed in Tab. 1. The operators 
in Tab. 1 define the general effective theory of dark matter-nucleon interactions mediated by 
a heavy spin-1 or spin-0 particle. The isoscalar and isovecotor coupling constants in (2.1), 
and cl respectively, are related to the coupling constants for protons (c^) and neutrons (c^) 
as follows: c| = (c° + c^)/2, and = (c° — c],)/2. These constants have dimension mass to 
the power —2. 

Assuming one-body dark matter-nucleon interactions, the most general Hamiltonian 
density for elastic dark matter-nucleus interactions is given by 

«T(r) = ^ (2.2) 

i=l T=0,1 k 


which is the sum of A identical terms. The index i in Eq. (2.2) identifies the zth-nucleon in 
the target nucleus, A is the mass number, and r is the relative distance from the nucleus 
centre of mass to the dark matter particle. 

The Hamiltonian density T-Ltj depends on the nuclear charges and currents, and admits 
the following coordinate space representation [16] 
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^Following [20], we neglect the operator O 2 = (v^)^, since it is quadratic in v^, and it cannot be a leading 
operator in effective theories. 
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where ai/2 and r* represent the the spin operator, and the position in the nucleus centre of 
mass frame of the ith-nucleon in the target nucleus, respectively. In Eq. (2.3) we use the 
definitions 
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where v;)^ = v-*- — v]^, and irj^ is an operator acting on the ith-nucleon space coordinate [16]. 
Explicit coordinate space representations for the operators q, v]^, and v;)^ can be found 
in [ 12 ]. 

Eq. (2.3) accounts for all possible ways of coupling dark matter to the nuclear charges 
and currents. Dark matter couples to the nuclear vector and axial charges through the 
operators Iq and and it couples to the nuclear spin, convection and spin-velocity cur¬ 
rents through the operators I 5 , \\j and 1^, respectively. Notice however, that the nuclear 
axial charge does not contribute to scattering cross-sections when nuclear ground states are 
eigenstates of P and CP. 

The operators in Tab. 1 contribute to Eq. (2.4) through specific combinations of q and 
v;)^. Notably, the operators Oi and O 4 generate the only momentum and velocity independent 
terms in fj, 1 ^, and 1 ^. 


2.2 Dark matter scattering from nuclei 

We use the Hamiltonian density (2.3) to calculate the amplitude Mnr for dark matter 
scattering from nuclei in the Sun. Summing (averaging) \Mnr\‘^ over final (initial) spin 
configurations one finds [16] 
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where the index k extends over (pairs of) nuclear response operators defined below in 
Eq. (2.8). 

The 8 dark matter response functions depend on matrix elements between eigen¬ 
states of of the operators (2.4). They are functions of q'^/rrij^ and = w'^ — (/^/(d/x^), 
where ixt and w are the reduced dark matter-nucleus mass and relative velocity, respectively. 
We list the functions in Appendix B. 

The 8 isotope-dependent nuclear response functions W'^'^ in Eq. (2.5) depend on the 
nuclear matrix elements of the charge and current operators in Eq. (2.3). For a given pair of 
nuclear operators and the nuclear response function , k = AB, is defined 

as follows 

WXW) = Y.{J,T,Mt\\ Al-Aq) ||J,r,MT)(J,T,Mr|| B^.^^iq) ||J,T,Mr), ( 2 . 6 ) 

L 

where the ket | J, T, Mt) represents a nuclear state of spin J, isospin T, and isospin magnetic 
quantum number Mt, and Mj is the nuclear spin magnetic quantum number. Here we use 
the notation WYb = , for A = H. The reduction operation in Eq. (2.6) is done via the 

Wigner-Eckart theorem: 

{J,Mj\ALM;rYMj) = ^ {J\\AtAU) ■ (2.7) 

In Eq. (2.5), the operators Alm;t and Blm-,t correspond to the following nuclear response 
operators 
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where M^(g'ri) = jL(gri)Y^^(HrJ , and MiMiq^i) = jL{qri)YLMAri)- The vector spher¬ 
ical harmonics, YAnAn), are defined in terms of Clebsch-Gordan coefficients and scalar 
spherical harmonics: 


= Y,{L'ml\\L'lLM)YL'A^Ae\ , ( 2 . 9 ) 
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where e;^ is a spherical unit vector basis. The 6 nuclear response operators in Eqs. (2.8) arise 
from the multipole expansion of the nuclear charges and currents in Eq. (2.3). The multipole 
index L in Eq. (2.8) is bounded from above: L < 2J. 

The isotope-dependent nuclear response functions must be calculated through 

detailed nuclear structure calculations for the most abundant elements in the Sun. In our 
analysis of dark matter annihilation signals from the Sun, we adopt the nuclear response 
functions derived in [12] through numerical shell model calculations carried out for the 16 
most abundant elements in the Sun. 

Combining Eqs. (2.5), (2.6) and (B.l), we can finally write the differential cross-section 
for dark matter scattering from nuclei of type i and mass rrii 
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which is in general a function of the nuclear recoil energy Er = g^/(2mj), and of the dark 
matter-nucleus relative velocity. 


3 Dark matter annihilation signals from the Sun 

Given the scattering cross-section (2.10), we now review how to calculate the rate of dark 
matter capture by the Sun, and the neutrino-induced muon flux at neutrino telescopes from 
dark matter annihilation in the Sun. 


3.1 Capture 


Weakly interacting dark matter particles of the Milky Way dark matter halo have a finite 
probability of elastically scattering from nuclei in the Sun. A dark matter particle becomes 
gravitationally bound to the Sun, when traveling through the solar medium it scatters to a 
velocity smaller than the local escape velocity. The rate of scattering from a velocity w to 
a velocity less than the escape velocity v{R) at a distance R from the Sun’s centre is given 
by [3] 


(w) = ^ Tiiw 0 



^ J J Eu^/w^ dER 


(3.1) 


where E = m^w‘^/2, is the dark matter particle initial kinetic energy, dui/dER is the dif¬ 
ferential cross-section for dark matter scattering from nuclei of mass and density in the 
Sun ni{R), and u is the velocity of the dark matter particle at i? —)■ oo, where the Sun’s 
gravitational potential is negligible. depends on the radial coordinate R, in that 

w = + v{R)‘^. The sum in the scattering rate (3.1) extends over the most abundant 

elements in the Sun, and the dimensionless parameters and are defined as follows 
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For a population of galactic dark matter particles with speed distribution at infinity given 
by f{u), the differential rate of capture per unit volume is given by [3] 
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Integrating (3.3) over a sphere of radius Rq, where Rq is the solar radius, one obtains the 
total rate of dark matter capture by the Sun 
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In this work, we consider the 16 most abundant elements in the Sun, i.e. H, ^He, ^He, 

14n, 16q^ ^"^Mg, ^’^Al, ^®Si, ^^S, ^‘^Ar, '^^Ca, ^®Fe, and ^®Ni, using the densities 

ni{R), and the velocity v{R) implemented in the darksusy code [39]. At the same time, we 
assume the so-called standard halo model [40], with a Maxwell-Boltzmann speed distribution 
for /(n), and a local standard of rest velocity of 220 km s“^ [41-43]. 

3.2 Annihilation 

Once captured, dark matter particles undergo subsequent scatters and sink into the centre 
of the Sun. In a variety of models [4-7], in the Sun dark matter annihilates into Standard 
Model final states, producing a flux of energetic neutrinos escaping the solar medium, and 
potentially observable at neutrino telescopes. 

The differential neutrino flux from dark matter annihilation in the Sun is given by [4] 


dK ~ 47rZl2 X 


(3.5) 


/ 


where E^, is the neutrino energy, the branching ratio for dark matter pair annihilation 
into the final state /, D the distance of the observer to the center of the Sun, and dNl/dE^^ 
is the neutrino energy spectrum at the detector from dark matter annihilation into the final 
state /. 

The annihilation rate, T^, in Eq. (3.5) is defined as follows 


Ta = - / d^xn^(x) (fJann^^rel) , 


(3.6) 


where n(x) is the dark matter space density in the Sun, and (uanni’rei) the thermally averaged 
dark matter annihilation cross-section times relative velocity. Eq. (3.6) implies the relation 
Ta = C'aA^/ 2, where N is the number of dark matter particles in the Sun. The constant 
Ca is given by 
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and no equal to the dark matter density at the core of the Sun. 

The number of dark matter particles N obeys the following differential equation in the 
time variable t 

N = C- CaN^ , (3.9) 


which admits the general solution 
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Element 

Average mass fraction 

Element 

Average mass fraction 

H 

0.684 

^4Mg 

7.30x10-4 

^He 

0.298 

"Ul 

6.38x10-5 

^He 

3.75x10-4 

^«Si 

7.95x10-4 

12g 

2.53x10“^ 

32S 

5.48x10-4 

I^N 

1.56x10“^ 

4UAr 

8.04x10-5 

^0 

8.50x10“^ 


7.33x10-5 


1.92x10“^ 

5B|r 

1.42x10-5 


3.94x10“^ 


8.40x10-5 


Table 2. List of average mass fractions for the 16 most abundant elements in the Sun as implemented 
in the darksusy program [39] (see also [44]). 


In our analysis we assume -v/CCU^o S> 1, where is the age of the Sun. This assumption 
implies T^ = (7/2, which corresponds to equilibrium between capture and annihilation of 
dark matter in the Sun. 

Neutrinos from dark matter annihilation in the Sun can be detected at neutrino tele¬ 
scopes observing an upward muon flux induced by charged-current neutrino interactions with 
nuclei in the material surrounding the detector. The differential neutrino-induced muon flux 
at the detector is given by 


dE^ 



dE'^V{E^,E'^-X) 


dacc{Eu,E'^) 


(3.11) 


where E^^ is the experimental threshold energy, Nt is the number of nucleons per cubic 
centimeter, A is the muon range, 'P{E^, E'^] A) is the probability for a muon of initial energy 
E'^ to have a final energy E^ after traveling a distance A inside the detector, and dacc/dE'^ 
is the weak differential cross-section for production of a muon of energy E'^. 

We evaluate Eq. (3.11) using muon yields generated by WimpSim [11], and pre-tabulated 
in darksusy [39]. Notice that WimpSim also accounts for the angular dependence of the upward 
muon flux, not included in Eq. (3.11) for simplicity. At the same time, WimpSim assumes a 
Gaussian space distribution for dark matter in the Sun, which is not generically true when 
the dark matter-nucleon scattering cross-section is momentum or velocity dependent [36]. 
Departures from this assumption are however negligible for a dark matter particle mass 
larger than about 30 GeV (see Eig. 3 in [36]). 

In computing the neutrino-induced muon flux, we calculate the rate of dark matter 
capture by the Sun, and the relevant dark matter-nucleus scattering cross-sections using our 
own routines and nuclear response functions [12]. We then compare our predictions with 
neutrino telescope observations as explained in the next section. 


4 Model independent analysis 

In this section we compare the general effective theory of one-body dark matter-nucleon 
interactions mediated by a heavy spin-1 or spin-0 particle with current neutrino telescope 
observations. We focus on data collected by the 79-string configuration of the IceGube/Deep- 
Gore observatory [9], and during 3109.6 days of SUPER-K dark matter searches [8]. Both 
experiments report limits on the neutrino-induced muon flux from dark matter annihilation 


- 8 - 




















Figure 1. Nuclear response functions WU and W^T, as a function of the momentum transfer q, for 
T ^ , and for the 16 most abundant elements in the Sun, as well as for Xe and 1. Conventions for 

colors and lines are those in the legends. 


above a given energy threshold, and below a given angle from the direction of the centre of 
the Sun. In our analysis, we use Tab. 1 from [9] as IceCube/DeepCore data, and Tabs. 1 
and 2 from [8] as SUPER-K data. 

We use darksusy to compute the muon yields at the detector, and our routines and 
nuclear response functions (see Figs. 1, 2 and 3) to calculate the rate of dark matter capture 
by the Sun. For the dark matter primary annihilation channel, we consider two extreme 
scenarios corresponding to dark matter pair annihilation into W~^W~, as for models with 
hard annihilation spectra, and to dark matter pair annihilation into bb, as for models with 
soft annihilation spectra. 

For rriy. < 80.3 GeV, the W~^W~ channel is not kinematically allowed. Following [8, 9], 
in our IceCube/DeepCore hard spectrum analysis we replace the W~^W~ channel with the 
hnal state t~^t~. We instead neglect the W^W~ channel in the hard spectrum analysis of 
the SUPER-K data. 

Demanding that the predicted neutrino-induced muon flux does not exceed its experi¬ 
mental 90% confidence level upper limit, we derive exclusion limits at the 90% confidence level 
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on the isoscalar and isovector coupling constants corresponding to the interaction operators 
in Tab. 1. 

For each interaction operator in Tab. 1, we consider the corresponding isoscalar and 
isovector coupling constants separately. In the figures, we report our 90% confidence level 
exclusion limits as a function of the dark matter particle mass, varying in the range 10 
GeV - 10 TeV. 

We also compare our limits from neutrino telescope observations with the 2D 90% cred¬ 
ible regions that we obtain from a Bayesian analysis of the LUX and COUPP direct detection 
experiments. The details of our COUPP and LUX analysis are discussed in Appendix A. 

4.1 Momentum/velocity independent operators 

We start with an analysis of the momentum and velocity independent interaction operators 
0\ and O 4 . The operator 0\ generates a constant spin-independent cross-section < 7 ®^ 
for dark matter interactions with protons (neutrons) given by 

SI ^ ^ . 

P TT 4 ’ 
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CY — C 
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Figure 3. Same as for Fig. 1, but now for the nuclear response functions and . 

Analogously, the operator O 4 generates a constant spin-dependent cross-section 
for dark matter-proton (dark matter-neutron) interactions given by 

_SD _ f^NJxOx + 1 ) \4 + ciP . _SD _ l^%jxUx + 1 ) |C 4 “ 

- Att 4 ’ “ dTT 4 • 

In the expressions above, fijv = my-mN/{m^ + mN) is the reduced dark matter-nucleon mass, 
and is the dark matter particle spin. For definitiveness, in the calculations we assume 
Jx = 1 / 2 . 

Fig. 4 shows the 90% confidence level exclusion limits that we obtain on the coupling 
constants c)* (top left panel), c\ (top right panel), C4 (bottom left panel), and C4 (bottom 
right panel). Solid red (blue) lines correspond to the hard annihilation spectrum analysis of 
the IceCube/DeepCore (SUPER-K) data, whereas dashed red (blue) lines represent our soft 
annihilation spectrum analysis of the IceCube/DeepCore (SUPER-K) data. Eor comparison, 
we also report 2D 90% credible regions obtained from LUX and COUPP as explained in 
Appendix A. 

Our exclusion limits on the coupling constants and C4 (top left and bottom left panels 
in Pig. 4) agree with those in [8, 9] at the few % level for small values of m^. Eor rriy. ~ 10 TeV, 


- 11 - 



































































[GeV] 



[GeV] 



[GeV] 



[GeV] 


Figure 4. Exclusion limits on the isoscalar and isovector coupling constants corresponding to the 
operators Oi and O 4 . Limits are presented at the 90% confidence level. Solid red (blue) contours 
correspond to an analysis of the IceCube (SUPER-K) data which assumes dark matter pair annihi¬ 
lation into W~^W~ (or into r+T“, see text at the beginning of Sec. 4 for more details). Dashed red 
(blue) contours refer to an analysis of the IceCube (SUPER-K) data which assumes dark matter pair 
annihilation into bb. For comparison, we also report the 2D 90% credible regions that we obtain from 
LUX (yellow) and COUPP (green). Coupling constants are expressed in units of my = 246.2 GeV. 


our findings differ from those of [8, 9] by up to 1 order of magnitude in the cross-sections, 
as a consequence of using different nuclear response functions. Here we adopt the nuclear 
response functions recently derived in [12] through numerical nuclear structure calculations. 
In contrast, in previous analyses approximate exponential form factors have often been used. 


- 12 - 












[GeV] 



[GeV] 



[GeV] 



[GeV] 


Figure 5. Same as for Fig. 4, but now for the operators O 3 and O 5 . 


4.2 Momentum/velocity dependent operators 

We now move on to describe our analysis of the momentum/velocity dependent dark matter- 
nucleon interactions operators. All interaction operators in Tab. 1 but Oi and O 4 depend 
on q, on v-*-, or on both. 

For momentum or velocity dependent operators exclusion limits are weaker than for 
Oi and O4. The reason is that for operators Ok, /c 7 ^ 1,4, the rate of dark matter capture 
by the Sun is suppressed by powers of and = w'^ — in Eq. (B.l). For 

instance, the operator Oq generates a term proportional to in , and the operator On 
generates a term proportional to q^ in . Nevertheless, dark matter can scatter in the 
Sun with a larger relative velocity w (and momentum transfer) than in terrestrial detectors, 
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Figure 6. Same as for Fig. 4, but now for the operators Oq and O 7 . 

where the scattering velocity is less then about 800 km s“^. Indeed, being + v^{R) 

and 620 km s“^ < v{R) < 1400 km s“^ [39], then w > 620 km s“^. Therefore, since the 
scattering cross-section (2.10) is evaluated at w, and not at u, velocity dependent scattering 
cross-sections weaken the exclusion limits from neutrino telescopes less dramatically than 
those from direct detection experiments. 

For momentum or velocity dependent operators, exclusion limits also depend on nuclear 
response functions different from WJJ , and W^, . For reference, in Figs. 1, 2, and 3 

we plot W]j', wit!, W^r\ and 1F|7', with r = r', for the 16 isotopes in 

Tab. 2. We do not report here and r ^ r' components of the other 

nuclear response functions, since they are not relevant when the interference between different 
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Figure 7. Same as for Fig. 4, but now for the operators Os and Og. 

operators, or between isoscalar and isovector components of the same operator is neglected, 
as in the present analysis. In [12] we list all nuclear response functions in Eq. (2.6) for the 
16 elements in Tab. 2 in analytic form. 

In general, exclusion limits depend on the solar composition. For reference, in Tab. 2 we 
report the average mass fractions of the 16 most abundant elements in the Sun as implemented 
in darksusy. 

Figs. 5, 6, 7, 8, 9, and 10 show the exclusion limits on the isoscalar and isovector 
coupling constants corresponding to the operators Ok, k = 3, 5, 6 ..., 15 that we find from 
IceCube/DeepCore and SUPER-K data. These figures reveal that the operators most severely 
constrained by current neutrino telescope observations are On, On, Og and O 3 (besides the 
operators Oi and O 4 discussed in the previous section). Notably, our current exclusion 
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Figure 8 . Same as for Fig. 4, but now for the operators Oio and On. 

limits on the interaction On = iS^ ■ (q/mjy), are more stringent that those set on the 
familiar O 4 = ■ Sjv interaction. We focus here on isoscalar interactions only. This result 

is expected [12], in that On is independent of the transverse relative velocity operator, and 
of the nucleon spin, which implies a large capture rate proportional to A^. 

4.3 Comparison with LUX and COUPP 

We now compare our limits on the constants from neutrino telescopes with the limits that 
we obtain from LUX and COUPP. 

For dark matter-nucleon interaction operators mainly coupling to the nuclear response 
operators and dark matter direct detection experiments tend to set stronger 

limits on the coupling constants c^. The reason is that and with some restriction 
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Figure 9. Same as for Fig. 4, but now for the operators O 12 and O 13 . 

also favor heavy elements. Indeed, in the g —)■ 0 limit, measures the nuclear 

mass number, and the content of nucleon spin-orbit coupling in the nucleus, which 

is large for heavy elements with orbits of large angular momentum not fully occupied [16] . 
Elements in the Sun are hence disfavored compared to Xe or I, which respectively compose 
the LUX and (partially) COUPP experiments. Interaction operators of this type are Oi, O 3 , 
O 5 , Og, On, O 12 and O 15 . 

The operators O 5 and Og also couple to the nuclear response operator which 

in the low momentum transfer limit measures the nucleon angular momentum content of 
the nucleus. The nuclear response operator Alm;t is generated by the nuclear convection 
current in Eq. (2.3). It favors nuclei with an unpaired nucleon in a non s-shell orbit, like for 
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Figure 10. Same as for Fig. 4, but now for the operators Ou and O 15 . 


instance and in the Sun’s interior. The coupling of the operators O 5 and Os to the 
nuclear response operator Alm;t strengthens our limits on C 5 and Cg from neutrino telescope 
observations. 

The remaining operators, namely O 4 , Oq, O 7 , Og, Oig, O 13 , and O 14 mainly couple 
to the nuclear response operators and which in the low momentum transfer 

limit measure the nucleon spin content of the nucleus. Within this group of interaction 
operators, of particular interests are the operators Oq and O 7 . For the operators Oq and 
O 7 , we find that neutrino telescopes can set stronger limits on Cg and C 7 than dark matter 
direct detection. Notably, for the operator Oj our exclusion limits on the coupling constant 
C 7 are stronger than those from COUPP and LUX in a large m-^ range, even within the 
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Figure 11. Isotope-dependent exclusion limits on the coupling constants cj (top left panel), c® (top 
right panel), Cy (bottom left panel), and Cy (bottom right panel) from a hard spectrum analysis of the 
IceCube/Deepcore data. Colored lines assume dark matter scattering from the specific isotope in the 
legend. Black lines correspond to total exclusion limits, i.e. including all isotopes. For comparison, 
in each panel we also report the 2D credible region that we obtain from the LUX experiment (yellow 
thick line). Limits are presented at the 90% confidence level. Though spin-independent, we also 
include the isoscalar component of Oi for reference. We also include the isovector component of O 7 , 
since limits from neutrino telescopes on this operator are particularly competitive. 


soft annihilation spectrum analysis. The competitive limit on from neutrino telescopes 
depends on the term in , which is less penalizing for neutrino telescopes than for 
direct detection experiments. In addition, it also depends on the abundance of Hydrogen 
in the Sun, which for O 7 determines rate of dark matter capture, and on the large value of 
W£T for this element. 


4.4 Spin-dependent operators and the heavy elements in the Sun 

We conclude this section exploring the spin-dependent dark matter-nucleon interaction op¬ 
erators more in depth. In particular we show that only for the spin-dependent operators 
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Figure 12. Same as for Fig. 11, but now for the coupling constant Cg, Cg, Cg, and c°g. 


O 4 and Oi Hydrogen is the most important element in the exclusion limit calculation. For 
all other spin-dependent operators, i.e. O3, Oq, Og, Oig, Oi2; O13, and O15, heavier 
elements are significantly more important. This conclusion is illustrated in Figs. 11, 12, and 
13, where we report exclusion limits obtained assuming that dark matter scatters from single 
elements in the Sun. Conventions for colors and lines are those in the legends. For simplicity, 
we focus here on the isoscalar components only. 

For the operators O 3 , O 12 and O 15 the most important element is ®®Fe, since they 
mainly couple to the nuclear response operator which favors elements with large A 

(and in the Sun ^®Fe is more abundant than ®®Ni). 

For the operators Oq, Og, Oiq, O13 and O 14 , exclusion limits mostly depend on 
since in this case dark matter scattering from Hydrogen is suppressed for > 20 GeV. 
Indeed, for rrii equal to the Hydrogen mass, the integral in Eq. (3.1) is largely dominated by 
small recoil energies, where the differential cross-section for these operators is particularly 
small. Notice also that the operators Oq, Og, Oiq, O13 and O14 mainly couple to the nuclear 
response operators that is characterized by fairly large W^T and 
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Figure 13. Same as for Fig. 11, but now for the coupling constant 0 ^ 2 , cjg, and 0 ^ 5 . 


W^J,' nuclear response functions. 

5 Conclusions 

We derived exclusion limits on the coupling constants of the general effective theory of one- 
body dark matter-nucleon interactions using current IceCube/DeepCore and SUPER-K ob¬ 
servations. In this study, we exploited new nuclear response functions computed in [12] 
through nuclear structure calculations for all dark matter-nucleon isoscalar and isovector in¬ 
teraction operators in Tab. 1, and for the 16 most abundant elements in the Sun. Exclusion 
limits were presented at the 90% confidence level, and separately assuming dark matter pair 
annihilation into W^W~ and bb. 

We found that the most severely constrained interaction operators are On = iSy- ■ 
(q/m 7 v), O12 = • (Stv X V-*-), Og = • v-*- and O 3 = iS^ ■ [i^/'mN) x v-*-], besides 

the familiar spin-independent and spin-dependent operators Oi = and O 4 = • Sat, 
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respectively. For each operator we found a physical interpretation for the relative strength 
of the corresponding limit. 

We compared our limits on the coupling constants from neutrino telescopes with bench¬ 
mark dark matter direct detection exclusion limits. We found that the operator Of = S^r • v*- 
is significantly more constrained by neutrino telescopes than by LUX and COUPP. This 
conclusion is mainly related to the term in R^, , which is less penalizing for neutrino 
telescopes than for direct detection experiments, to the abundance of Hydrogen in the Sun, 
and to the corresponding large value of W^T . 

Another important result of this work is to show that Hydrogen is not the most im¬ 
portant element in the exclusion limit calculation for the majority of the spin-dependent 
operators in Tab. 1. As a consequence, nuclear structure calculations as those initiated 
in [12] appear to be key tools for model independent analyses of dark matter signals at 
neutrino telescopes. 
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A LUX and COUPP analysis 


We compare the general effective theory of one-body dark matter-nucleon interactions me¬ 
diated by a heavy spin-1 or spin-0 particle with the LUX and COUPP experiments in a 
Bayesian analysis based on the Likelihood function 


- ln£(d|mx, c) ~ c) + {2 - k) In c) + As 


— In 


{k‘^ - k) 


1 2 


(^B+ + jlB - 


+ constant. 


(A.l) 


In Eq. (A.l) k is the number of observed recoils in a given dataset d, c) is the number 

of predicted scattering events at a given mass m^, and coupling constant c, and As is the 
corresponding number of expected background events (with error us). The constant term in 
Eq. (A.l) corresponds to the arbitrary normalization of C, which we fix by imposing InT = 0, 
for As = k. 

Eor the COUPP Likelihood function we assume k = 2, k = 3 and k = 8 bubble nucle- 
ations above a threshold energy Eth of 7.8 keVnr, U keVnr and 15.5 keVnr, respectively [45]. 
For the 3 threshold configurations the expected number of background events is, respectively, 
/is = 0.8, /is = 0.7 and /is = 3. In the analysis we use an energy dependent experimen¬ 
tal exposure e characterized by: e(7.8 keV) = 55.8 kg-days, e(ll keV) = 70 kg-days and 
e(15.5 keV) = 311.7 kg-days [45]. For COUPP, we compute /ig as follows 

E f°° dlZT 

/ dERVTiEji, Eth)-7^ (A.2) 

T=C,F,1 
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where (HZt/^Er is the differential rate of scattering from the target material T. The prob¬ 
ability VTiEu, Eth) that an energy Er nucleates a bubble above -Eth is given by [45]: 


Vt{Er, Eth) = 1 - exp 


— OiT 


Er — Eth 

-Eth 


(A.3) 


We assume perfect efficiency for bubble nucleation for Iodine (oj —>■ +oo), we neglect scat¬ 
tering from Carbon (etc = 0), and marginalize over etp = OCOUPP) assuming for the latter a 
log-prior with mean 0.15. 

We construct the LUX Likelihood function as in [25], assuming k = 1 and jls = 0.64 ± 
0.16. We consider a Gaussian resolution for photoelectron detection with standard deviation 
<7PMT = 0.37, an exposure of 250x85.3 kg-days, and the experimental efficiency in Fig. 1 of 
Ref. [46] times 0.5, in order to account for the 50% nuclear recoil acceptance quoted by the 
LUX collaboration. 

Results are presented in terms of x% credible regions, i.e. portions of the parameter 
space containing x% of the total posterior probability, and where the posterior probability 
density function 7^(0|d) oc£(d| 0 ) 7 r( 0 )at any point 0 in parameter space inside the credible 
region is larger than at any external point. We assume a uniform prior probability density 
function 7 r( 0 ) for the logarithm of and of c. Prior and volume effects in dark matter 
direct detection are discussed in [25, 47-49]. 


B Dark matter response functions 


In this appendix we list the dark matter response functions that appear in Eq. (2.10). The 
notation is the same used in the body of the paper. 
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